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Abstract. We present part of the results from a spectro- 
scopic monitoring programme on a sample of AGNs, rela- 
tive to Narrow Line Seyfert 1 Galaxies: following the idea 
that Balmer-line variability can help discriminate among 
the possible models for these objects, and on the basis of a 
first monitoring project with long time-scales, we now add 
to the previous database new variability results, both on 
short and long time-scales. The comparison with a similar 
data-set of normal Seyfert l's suggests now a statistically 
different behaviour of NLSls with regard to the 1-year 
variability properties, and a clearer trend of NLSls to- 
wards weaker variability is suggested by the comparison 
between the 1-month variations of NLSls and those of 
the typical Seyfert 1 NGC 5548. Although these results 
do not take into account the possible role of the relative 
contribution of narrow components to the line fluxes, and 
of the continuum variability amplitudes (both of which 
can in principle be different in the two classes of objects), 
they are consistent with the model of a BLR which is on 
average larger in NLSls. 

Key words: Galaxies: active - Galaxies: nuclei - Galax- 
ies: Seyfert 



1. Introduction 

The widespread property of variability of the contin- 
uum and broad lines emitted by Active Galactic Nuclei 
has been used in recent years to obtain constraints on 
the structure and kinematics of the Broad Line Region, 
through intensive spectral monitoring and the use of rever- 



beration mapp ing techniques (Blandford & McK ee 1982, 
Peterson |l993| , Robinson |l994| , Alloin et al. |l995|) : the in- 
ferred sizes of the BLR range from a few light days to 
about a light month, and therefore direct information on 
the region through imaging techniques is not available. 

Send offprint requests to: M. E. Giannuzzo 



Narrow Line Seyfert 1 Galaxies (Osterbrock & Pogge 



1985 ) are a sub-class of Seyfert l's located at the lower 
end of the broad line width distribution: the FWHM of 
the permitted lines does not exceed 1000-1500 km s _1 in 
these objects. In most other properties NLSls are sim- 
ilar to 'normal' Seyfert l's: they have different widths 
of permitted and forbidden lines, intense Fen lines and 
high ionization transitions, normal line ratios and average 
luminosities. They present, on average, lower equivalent 



widths compared to Seyfert l's (Osterbrock & Pogge 1985, 
Goodrich 1989), but again this is the extension to low 



FWHM of a trend observed throughout the whole Seyfert 
1 population. NLSls are efficiently found in soft X-ray se- 
lected Seyfert 1 samples, in which they represent ~ 16-50% 



of all objects (Stephens 1989 , Puchnarewicz et al. 1992 ), 
compared to the ~ 10% found in optically-selected sam- 
ples. Boiler et al. ( 1996j ) found that NLSls have generally 
much steeper soft X-ray (ROSAT) continuum slopes than 
normal Seyfert l's, and present in a few cases rapid soft X- 
ray variability. The 2-10 keV spectra have also been found 
to be steep er in NLSls than in other Seyfert l's (Brandt 
et al. |1997| ). 

The line-of-sight component of the velocity of the 
broad line emitting gas may be unusually small for several 
reasons, among which a lower mass of the central black 
hole, projection effects due to gas motions in a plane (e.g. 
that of a disk) observed almost pole-on, a larger BLR size, 
a lack or obscuration of the inner, high velocity regions of 
the BLR. The explanation of this behaviour can therefore 
provide insight on the more general problem posed by the 
great diversity among broad emission line profiles. 

Permitted line variability, being directly related to the 
size, geometry and kinematics of the Broad Line Region, 
is potentially very useful to discriminate among some of 
the mentioned hypotheses: a lack of widespread variability 
in NLSls would in fact suggest that the visible BLR is lo- 
cated at a higher distance from the central engine than in 
normal Seyfert l's; similar variability properties in NLSls 
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and Seyfert 1 's could instead imply a smaller central mass 
or an anisotropic kinematic structure for the BLR. 

Since no NLS1 has ever been intensively monitored, 
and virtually no information was available on the pres- 
ence of variability in this class of objects when the project 
started, we performed a preliminary simple observational 
programme aimed at determining whether variability is a 
common characteristic of NLSls as a class (Giannuzzo & 
Stirpe 1995, Giannuzzo & Stirpc 1996 - hereafter Paper 
I): we adopted a statistical approach, by evaluating the 
optical line flux percentage variations on a time-scale of 
1 year for a sample of NLSls, and comparing our results 
with those of an existing similar data-set of normal Seyfert 
l's. We did not find evidence for a weaker or less common 
variability in Narrow Line Seyfert l's, and discussing the 
competing models we concluded that a larger BLR or a 
smaller black hole mass could be the most probable expla- 
nations; we stressed that a short-term monitoring program 
could help discriminate between the two possibilities. This 
kind of observations could also address the problem of 
the generality of reverberation mapping results, since until 
now only AGNs with high variability and normal line pro- 
files have been monitored, perhaps discriminating against 
different types of BLRs (Robinson 1995 ). 

In this work the observations reported in Paper I are 
added to new data on NLSls from a spectroscopic moni- 
toring programme, conducted at the Bologna Observatory 
on a larger sample of AGNs, which lasted a few years: this 
following step of the project aims therefore at collecting 
information on the variability properties of a greater num- 
ber of objects on different time-scales, though not in much 
detail. We have thus enlarged the previous data-set on 1- 
year percentage variations of the optical lines (this time 
both H/3 and Ha have been considered), allowing in this 
way a more accurate statistical comparison than before; 
furthermore, we evaluated the short-term variations of the 
line fluxes with the aim to compare them, in the case of 
Ha, with those of some monitored Seyfert l's, while in 
the case of H/3 we used as a comparison object the best 
studied Seyfert 1, NGC 5548, the light curve of which has 
been well sampled for several years (Korista et al 
and references therein). 

Section 2 describes the observations and data reduc- 
tion and analysis, section 3 presents our results and the 
statistical comparisons between NLS1 and Seyfert 1 vari- 
ability properties; finally, section 4 includes a discussion 
on the results and our conclusions. 



2. Observations, data reduction and analysis 

The sample of AGNs chosen for the monitoring campaign 
includes 5 objects classified (sometimes 'marginally') as 
Narrow Line Seyfert l's, 2 of which were included also 
in the sample observed at ESO (Paper I): more specifi- 
cally, the NLS1 group includes th e pro totype of the class, 
I Zwicky 1 (e.g. Halpern & Oke |l987| ), Akn 564 (Oster- 



1995 



brock & Shuder |l982j ), NAB 0205+024, Mkn 335 and PG 
1211+143 (classified by us on the basis of published spec- 
tra according to the Osterbrock & Pogge 1985 criteria). 
The last two objects are considered as 'marginal' because 
their permitted line FWHM is slightly larger than the ~ 
1000 km s" 1 limit, being 1200 km s _1 in Mkn 335 (unpub- 
lished data collected by the authors) and 1600 km s _1 in 
PG 1211+143 (Stirpe |1990| ), but in both cases the narrow 
line had been subtracted, which increased the measured 
line FWHM; furthermore, they share other properties with 
typical NLSls, such as intense Fen emission and weak 
forbidden lines. A sample continuum-subtracted spectrum 
for each of the observed NLSls is plotted in Figure [j]. 

The observations were carried out at the Cassini 1.52 m 
telescope operated by the Bologna Observatory, using two 
different spectrographs. A log of the observations is given 
in Table |l|. In the first part of the monitoring campaign 
(October 1991 - March 1994) the spectra were obtained 
with a Thomson 7882 CCD (576x384 pixels) mounted on 
a Boiler & Chivens spectrograph. The 350 1/mm grating 
used at first order gives a coverage of 2500 A with a reso- 
lution of 4.5 A/pixel. The grating angle was adjusted for 
each object in order to cover, when possible, both the Ha 
and H/3 spectral regions. In the summer of 1994, a new all- 
transmitting spectrograph named BFOSC (Bologna Faint 
Object and Spectrograph Camera) replaced the Boiler & 
Chivens. The new instrument, working at the Cassegrain 
focus of the telescope as a focal reducer, is based on a 
EFOSC concept (Enard & Delabre |1982|). The observa- 



tions after July 1994 were therefore carried out with a 
Thomson coated CCD (1024x1024) mounted on BFOSC, 
with a projected pixel size of 0.56 arcsec. The used grism 
gives a larger but fixed spectral range (4000-7850 A) and 
a resolution of 4 A/pixel. 

The slit width was 2.5 arcsec, matching the typical 
seeing of the site, with the exception of the March 1995 
observation of NGC 5548, for which a double width was 
chosen due to bad seeing conditions, and a few observa- 
tions of Akn 564 and Mkn 335, for which the good seeing 
ensured that also with a slit 2.0 arcsecs wide all the flux 
entered the instrument. 

The observations were performed with standard pro- 
cedures, including the acquisition of bias and dome flat 
field frames, for subsequent background and pixel to pixel 
sensitivity correction, and of calibration lamp frames for 
wavelength calibration; standard stars were observed to 
perform flux calibration. 

Together with the data on the NLSls of our sample, 
we were interested in the NGC 5548 observations, to be 
added to literature data in the variability comparison as 
a typical Seyfert 1 galaxy. 

The data reduction was performed with standard 
IRAF tasks. As is usual in variability studies, to obtain 
line fluxes which are accurate within a few percent we 
also intercalibrated our spectra to correct for slit losses, 
by using the fact that the strong forbidden lines emitted 
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Galaxy a (1950) 5 (1950) my z Date of observation Julian date Integration time 

(JD - 2400000) (minutes) 



Mkn 335 00 03 45.2 +19 55 29 13.85 0.025 



I Zwicky 1 00 50 57.8 +12 25 20 14.03 0.061 



NAB 0205+024 02 05 14.5 +02 28 42 15.40 0.155 



PG 1211+143 12 11 44.8 +14 19 53 14.63 0.085 



NGC 5548 14 15 43.5 +25 22 01 13.73 0.017 



Akn 564 22 40 18.3 +29 27 48 14.16 0.025 
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Fig. 1. Representative continuum-subtracted spectra of the monitored Narrow Line Seyfert 1 objects. 



by the Narrow Line Region remain constant on very long 
time-scales, and can therefore be used as standard candles. 
Since the [S n] A6717.0, A6731.3 A lines are quite weak and 
not well resolved in all spectra, we used the strong [O III] 
A4958.9, A5006.8 A transitions for the calibration in the 
whole observed spectral region: as a result, the correction 
is much more accurate in the H/3 portion of the spectrum 
than in that of Ha, and in a way not easy to quantify. The 
correction parameters (wavelength shift and scale factor) 
to be applied to all the spectra of each object, with refer- 



ence to one chosen spectrum, are automatically calculated 
by a code (van Groningen & Wanders 1992) which makes 
use of a x 2 minimum research procedure on the difference 
spectrum in the wavelength range of the forbidden lines. 
We found that after the correction the integrated fluxes 
of the two [O III] lines for each object differ by less than 
10% from their average value. 

A possible problem of this method, however, is the 
fact that different seeing conditions may cause different 
portions of the NLR and the BLR to enter the slit in 
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observations at different epochs, therefore leading to errors 
in the calibration. In Paper I we performed several tests 
of the data against this source of error, concluding that, 
within the uncertainties intrinsic to the method, this kind 
of problem does not appreciably affect the reliability of 
the internal calibration. 

To the spectra thus obtained from the Loiano pro- 
gramme, we added (after having performed on them the 
usual internal calibration) two spectra of Akn 564 and 
Mkn 335 obtained at the ESO 1.52 m telescope in October 
1996 as part of another programme: the covered spectral 
range is in this case ~ 2400 A with a resolution of 4.6 
A/pixel, and the observations have been performed with 
a 2 arcsec slit. We also included some of the data from the 
ESO observations used in Paper I. Table |l| also includes 
information on these observations, thus giving a complete 
list of the data taken by the authors and used throughout 
this paper. 

3. Variability results and comparison with Seyfert 
l's 

The fluxes of the strongest optical lines (Ha and H/3) were 
measured from all our spectra, for each object and each 
observing date, by fitting a straight-line continuum under 
the lines and integrating the overlying fluxes. 

The uncertainties on these fluxes, and consequently 
on the calculated percentage variations, were evaluated as 
described in Paper I, the major contribution being the un- 
certainty on the internal calibration mentioned above; this 
contribution was then estimated by slightly perturbing 
the scale factors of the spectra and measuring the ranges 
within which no appreciable changes could be seen. The 
obtained percentage uncertainties on the line fluxes fall in 
the range 3-5%, and the errors on the flux variations we 
used for the following analysis are therefore approximately 
twice this amount. 

Figure || shows the total H/3 light curves of the ob- 
served NLSls: similarly to what we found for the 'ESO 
sample', most objects appreciably varied their line fluxes 
at least in some of the considered time intervals. The H/3 
light curve of NGC 5548, as given by our Loiano observa- 
tions only, is also plotted in the figure. 

In Table || the measured Ha and H/3 fluxes for the 
Loiano and ESO spectra are reported. Not having enough 
literature data on the Seyfert 1 short-term variability of 
other emission lines (such as Fen or Hen), we did not 
measure their fluxes and relative variations. We did not 
evaluate the continuum flux either, because it is expected 
to be contaminated by the host galaxy starlight in differ- 
ent amounts for different observations, due to the differ- 
ent seeing conditions, slit apertures and source distances; 
therefore it is not possible to determine the global contin- 
uum variability properties in a reliable way. 

The variability behaviour of individual Narrow Line 
Seyfert l's in the Loiano (+ ESO) sample is quantified 



in Table |[ which reports, for each object, the rms (er) 
of the H/3 flux throughout the whole campaign, and the 
fractional variation F var , defined (Clavel et al. 1991) as the 



ratio of the rms fluctuation to the mean line flux, corrected 
for the mean measuring error. Note that this last quantity 
is very similar (~ 5-9%) for all the NLSls, while it is 
much higher for NGC 5548, the only Seyfert 1 observed 
at the same telescope and approximately on the same time 
baseline, and the data of which are therefore also reported 
in the table for comparison. This is a first indication of a 
lesser degree of variability amplitude displayed by Narrow 
Line objects. 

To perform a direct quantitative comparison, we first 
evaluated, for each NLS1 in the 'Loiano sample', all the 
relative flux variations [] on time intervals ranging from ~ 
10 months to ~ 14 months; to these data we added the 
similar ones obtained from the previous ESO campaign 
(provided that there are no overlapping intervals, to avoid 
to include twice the same relative variation) and from 
the 1996 ESO observations mentioned above. We thus ob- 
tained a data-set of 28 points (i.e. relative flux variations) 
for H/3 and 27 points for Ha. 

As for the Seyfert l's, also in this case we enriched 
the previous data-set: first we added to the de Ruiter & 



Lub Se yfert 1 sample (de Ruiter & Lub |1986| , Lub & de 
Ruiter 1992; see also Paper I for details on the sample 



- line fluxes have been obtained in the same way as for 
ESO and Loiano spectra) the data of the typical Seyfert 1 
NGC 5548 taken from Loiano observations for the period 
since 1992. We then took some points from the Ha and 
H/3 NGC 5548 literature data that we collected to pe rform 



the short tim e-scale comparison (Wamsteker et al. 1990 



Korista et al. 1995 and references therein). The data-set 
includes in this case 135 points for H/3 and 88 for Ha. 

To compare the short-term variability properties of the 
two classes of AGNs, we chose a short average time inter- 
val, such that the relative variations measured on this scale 
would be both representative of the characteristics of the 
objects and numerous enough to allow as accurate a statis- 
tical comparison as possible. Since the Loiano light curves 
are not evenly sampled, and the time gaps are sometimes 
of the order of several months, we took all the relative vari- 
ations of the optical lines on time intervals ranging from 
~ 10 days to ~ 45 days, the average interval being around 
27 days. In this way we used all the data, maximizing the 
available information and obtaining finally 23 points for 
H/3 variations and 17 points for Ha ones, having added to 
the H/3 set also 2 variations calculated with ESO Paper I 
observations. 

A similar kind of information can be obtained, in the 
case of H/3, for the typical Seyfert 1 galaxy NGC 5548, 



1 Here and in the following we always refer to relative varia- 
tions taken from single measurements - only the 5 ESO points 
are the result of averaged spectra - and computed with respect 
to the flux mean value. 
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Table 2. Measured Balmer line fluxes (some points considered less reliable are indicated in parenthesis). 



Galaxy 



Mkn 335 



I Zwicky 1 



NAB 0205+024 
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1 Units are 10 15 erg s 1 cm 



Table 3. Global variability properties of the H/3 flux of each object within the Loiano + ESO monitoring campaigns 
(see text for details). 



Galaxy 



Classification a 1 {Fh/}) F„ ar (Fjj/j) 



Mkn 335 
I Zwicky 1 
NAB 0205+024 
PG 1211+143 
NGC 5548 
Akn 564 



NLS1 
NLS1 
NLS1 
NLS1 
Sy 1 
NLS1 



27 
21 

8 

40 
102 



8.9% 
5.4% 
6.6% 
6.5% 
24.6% 
8.3% 



1 Units are 10 15 erg s 1 cm 
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Fig. 2. Total H/3 light curves of the observed NLSls and of the Loiano observations of NGC 5548. 



which we take here as representative of the Seyfert 1 class 
as regards variability characteristics (although, as men- 
tioned in Section 1, the generality of its properties is to 
be determined): since in fact this object has been optically 
monitored for several years with high temporal resolution, 
it is possible to check its monthly variations by properly 
sampling its H[3 light curve for the period ~ 1989-1992 
(Korista et al. 1995 and references therein). The sampling 
has been performed by simply taking the observation clos- 
est in time to the epoch which is separated from the pre- 



vious point by a time interval equal to the average chosen 
one (in this way we also 'mimick' the uneven sampling of 
NLS1 data). To these observations, sampled so that the 
average time interval is of the order of 27 days, we also 
added the Loiano data of the same object which were not 
temporally overlapped to the literature light curve. We 
collected in this way a total of 56 points for H/3 varia- 
tions. In the case of Ha, we tried to collect literature data 



on several LAG (Robinson 1994) and AGN Watch mon- 
itored objects, again sampling the literature light curves 
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to obtain line flux relative variations on ~ 1 month time- 
scale. However, this last data-set resulted to be very in- 
homogeneus; on the other hand, the [Olll] correction of 
the Loiano data does not give for the Ha region a calibra- 
tion accurate enough to allow a reliable comparison, and 
therefore we did not use the Ha 1-month variations in the 
following analysis. 

Similarly to what we did in Paper I, a first statistical 
comparison between Seyfert l's and NLSls can be per- 
formed with these data by constructing, for each class 
of object, time-scale of variations, and emission line, a 
histogram including all the measured relative variations 
(in absolute value). We then compared the resulting his- 
tograms relative to Seyfert l's and NLSls for each group 
of parameters, as plotted in Figure ||a-c: the 1-year data 
obtained for Ha and H/3 variations seem now to show a 
trend of Narrow Line Seyfert l's towards weaker variabil- 
ity, and a clearer trend in this sense is suggested by the 
1-month H/3 results, where the NLS1 histograms display a 
steeper decline and go to zero at around 15-20% percent- 
age variability. 

A more quantitative way to evaluate the statistical be- 
haviour of the two classes of objects is the application of 
the Kolmogorov-Smirnov (K-S) test, which gives the prob- 
ability (P) that two continuous data samples are drawn 
from the same parent population, together with the maxi- 
mum distance (D) between the 2 cumulative probabilities; 
since P depends on the number of data points in each 
sample, larger samples have higher statistical significance, 
while more conservative results are obtained for smaller 
samples. This is the case, for example, of the H/3 line flux 
data reported in Paper I, for which the K-S test gives 
D ~ 0.2 and P ~ 58%: our conclusion that there was no 
evidence for a different behaviour was therefore justified 
at that time. 

The same test on the enlarged 1-year samples, however, 
gives now different results, i.e. D ~ 0.3 and P ~ 5.4% 
for H/3 and D ~ 0.3 and P ~ 3.9% for Ha; so it seems 
improbable that the 2 sets of data are taken from ob- 
jects with similar variability properties. A lower probabil- 
ity that their parent population is the same is obtained 
by the K-S test applied to H/3 1-month data: here, in fact, 
we got D ~ 0.4 and P ~ 0.8%. Note that neither the his- 
togram appearance nor the K-S test results for the short- 
term variations result substantially altered if we sample 
the NGC 5548 H/3 light curve in a slightly different way 
(e.g. starting from a different point), the probability P al- 
ways ranging between 0.5 and 0.9% and the distance D 
between 0.41 and 0.43. 

4. Discussion and conclusions 

As already mentioned in Section 1, we discussed the vari- 
ous interpretative models for Narrow Line Seyfert 1 spec- 
tra in Paper I, concluding that the possibilities of a smaller 
black hole mass or a BLR size which is larger on average 



are the most probable ones. The other interpretations are 
in fact in some way inconsistent with observational re- 
sults: for example, a projection effect explanation would 
imply too narrow a cone of visibility for NLSls to jus- 
tify their being more than 10% of the known Seyfert 1 
population; a depletion of the innermost part of the BLR 
would produce very low line equivalent widths, while the 
obscuration hypothesis is quite unlikely, since the results 
of X-ray observations all suggest that we have a direct 
view of the nucleus. 

The knowledge of the variability behaviour of these ob- 
jects, compared to other Seyferts, can therefore be crucial 
in the understanding of their nature, allowing us to esti- 
mate the BLR size. However, the analysis of variability 
data can be potentially influenced by some effects: first 
of all, if the relative contribution of the (constant) nar- 
row component to the total line flux is more important in 
NLSls than in normal Seyfert l's, then a global weaker 
variability is observed in the former objects, even if their 
broad component fluctuations are on average of the same 
amplitude as those of Seyfert l's. However, the evaluation 
of the ratio between narrow and broad component of Ha 
and H/3 in the Stirpe ( |1990| ) Seyfert sample shows that the 
relative contributions to the total line fluxes are not dra- 
matically different in the two classes of objects, the nar- 
row component being on average ~ 4% of the broad one 
in normal Seyfert l's and ~ 6% in NLSls; furthermore, it 
results to be 6% in NGC 5548, which assures that, at least 
in the short-term comparison with this object, our mea- 
sured relative variations are smoothed by the same factor. 
Another possible bias in the interpretation of the results is 
the unknown continuum variability which drives the line 
fluctuations: again, an intrinsically lower continuum vari- 
ability amplitude in NLSls would result in a weaker line 
variability, even if the BLR has a 'standard' size in these 
objects. Given the lack of other information (e.g. correla- 
tions between optical/UV variability amplitude and broad 
line width), we assume that the ionizing continuum has 
the same variability properties in these objects as in other 
Seyferts. In both cases, since a more detailed analysis of 
the data (involving the subtraction of narrow line com- 
ponents and host galaxy continuum from the spectra, to 
evaluate the continuum fluctuations) is beyond the scope 
of this paper, we will assume that the potential 'smooth- 
ing' effects are small compared to the one of which we 
found evidence. 

Based on the ensemble of data that we collected in this 
work, the weaker variability of NLSls, implied by the his- 
togram comparisons and the K-S test both for long- and 
short-term observations, seems to exclude that the broad 
line emitting gas is located as close to the centre as in nor- 
mal Seyfert l's, as assumed in the 'small black hole mass' 
model. It seems more probable, instead, that in these ob- 
jects the gas extends further out: its distance should be 
not so high as to completely smear out the ionizing contin- 
uum variations in the line reprocessing, but large enough 
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Fig. 3. Statistical comparison between the variability properties of Seyfert 1 galaxies (open histograms) and NLSls 
(shaded histograms): (a) H/3 flux variations over a time interval of ~ lyr; (b) same as (a) but for Ha variations; (c) 
H/3 variability on a monthly time-scale. 



to 'smooth' their amplitude more than normal BLRs do, 
producing in this way the statistically weaker variability 
that we observe. This is of course assuming that the phys- 
ical conditions (like the electron density or the column 
density) in the ionized gas emitting the broad lines are 
similar in NLSls and in normal Seyfert l's. 

The fact that we actually observed short-term varia- 
tions in these objects, though with small amplitude, im- 



plies that, in the 'extended BLR' picture, the emissivity 
radius of their BLR gas is probably closer to our estimate 
of r m i n (N LSI) ~ 3r m i n (Seyl), made in the case of this 
model in Paper I, than to the upper limits that can be set 
on the basis of the average time-scale of the observations: 
the 1-month histograms we constructed include in fact 
percentage variations measured on several time intervals, 
the mean of which is of the order of 27 days, but which also 
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include At ~ 10-20 days; the sample is not large enough 
to allow a more detailed evaluation of the variability at 
the level of a few days, without losing in statistical sig- 
nificance. Only through an intense monitoring campaign 
on one or more NLSls (one observation every few days) it 
would be possible to estimate the BLR size by measuring 
the lag between line and continuum light curves. 
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